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ABSTRACT 

We explore the role of complex multipolar magnetic fields in determining physical 
processes near the surface of rotation powered pulsars. We model the actual magnetic 
field as the sum of global dipolar and star-centered multipolar fields. In configurations 
involving axially symmetric and uniform multipolar fields, 'neutral points' and 'neu- 
tral lines' exist close to the stellar surface. Also, the curvature radii of magnetic field 
lines near the stellar surface can never be smaller than the stellar radius, even for very 
high order multipoles. Consequently, such configurations are unable to provide an effi- 
cient pair creation process above pulsar polar caps, necessary for plasma mechanisms 
of generation of pulsar radiation. In configurations involving axially symmetric and 
non- uniform multipoles, the periphery of the pulsar polar cap becomes fragmented 
into symmetrically distributed narrow sub-regions where curvature radii of complex 
magnetic field lines are less than the radius of the star. The pair production process 
is only possible just above these 'favourable' sub-regions. As a result, the pair plasma 
fiow is confined within narrow filaments regularly distributed around the margin of 
the open magnetic flux tube. Such a magnetic topology allows us to model the system 
of twenty isolated sub-beams observed in PSR B0943 + 10 by Deshpande & Rankin 
(1999, 2001). We suggest a physical mechanism for the generation of pulsar radio emis- 
sion in the ensemble of finite sub-beams, based on specific instabilities. We propose an 
explanation for the subpulse drift phenomenon observed in some long-period pulsars. 
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1 INTRODUCTION 

Observations and theory suggest that complex multipolar 
magnetic fields prevail near the surface of neutron stars, and 
play an important role in the physics of rotation powered 
pulsars. 

The complexity of surface magnetic fields is determined 
by the evolution of magnetic fields in neutron stars. Ac- 
cording to one scenario (e.g., Blandford, Applegate & Hern- 
quist 1983), the magnetic field is generated by currents 
flowing in the thin outer crust of the star, with thickness 
Ar <t. Rb, where Rs ~ 10^ cm is the neutron star radius. In 
this case, the superposition of magnetic multipoles of order 
I ~ Rs/Ar > 10 dominates the exterior field (Arons 1993). 
Assuming that the evolution of magnetic fields in isolated 
neutron stars is due to ohmic decay, Mitra, Konar & Bhat- 
tacharya (1999) show that the evolution of such high-order 
multipolar components {I < 25) is very similar to that of the 
dipolar field: the structure of the pulsar surface field does not 



vary much as the pulsar ages. An alternative approach by 
Ruderman (1991a,b,c) attributes the surface magnetic field 
evolution to changes in the magnetic field present in the 
core. Due to stresses, caused by the spinning down or the 
spinning up of the neutron star, strong enough to drive the 
crustal lattice beyond its yield strength, moving magnetized 
crustal platelets are created and account for 'sunspot-fike' 
clumps of magnetic field at the stellar surface. 

One first reason for introducing non-dipolar magnetic 
fields close to the stellar surface lies in the inability for the 
pure surface dipolar field to provide sufficient pair produc- 
tion, as required for the generation of pulsar emission. In ef- 
fect, in their classical paper, Ruderman & Sutherland (1975) 
state that the value of the radius of curvature of magnetic 
field lines, p, should be about Rs ~ lO'' cm near the stellar 
surface, in order to fulfill the conditions for a copious pair 
production process in the polar vacuum gap. However, such 
small curvature radii cannot be achieved assuming a pure 
dipolar magnetic field near the stellar surface: the contri- 
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bution of non-dipolar magnetic components, present in the 
close vicinity of the neutron star in the form of either mul- 
tipoles of order higher than the dipole one, or the 'sunspot- 
like' clumps, is necessary. These components are supposed 
to be strong enough to change the local topology of the 
magnetic field in the vicinity of pulsar polar caps. Arons 
& Scharlemann (1979) also point out that pair creation in 
pulsars with very long periods is impossible without invok- 
ing magnetic fields more complex than the dipolar one. In 
a thorough study of the space-charge limited flow of elec- 
trons in combined surface dipolar and quadrupolar magnetic 
fields, Barnard & Arons (1982) conclude that the curvature 
radii of resulting magnetic field lines can be much smaller 
than those of the pure dipolar field, and discuss the possible 
impact of this effect on both the pair creation process and 
the observations. 

Chen & Ruderman (1993) specify the conditions for a 
sufficient pair production to be triggered off, as required for 
the generation of pulsar emission. They conclude that only 
a very twisted (ps ~ Rs, where ps is the curvature radius 
at the stellar surface) and strong {Bb = 2 x 10^^ G) mag- 
netic field ensures that on the P — P diagram, the death 
line leaves all pulsars above itself, that is to say, 'permits' 
their existence. However, an exception is the recently dis- 
covered PSR J2144 - 3933 with period P = 8.5 s (Young, 
Manchester & Johnston 1999), which appears to lie even 
below this death line. With this latter observation in mind, 
Gil & Mitra (2001) study anew the pair creation process 
in a super-strong magnetic field. They conclude that all 
observed pulsars, including PSR J2144 — 3933, lie above 
the death lines deduced from their model, which adopts 
Bs/B^ ~ 100 (although still B^ < Sq) and ps ~ Q.lRs- 
Here, Sq = rr?(? jeh = 4.4 x 10^^ G is the critical magnetic 
field strength above which the photon splitting phenomenon 
may inhibit the pair formation process (Daugherty & Hard- 
ing 1983), while B^ and Bg are respectively the dipolar and 
the total - dipolar plus multipolar - magnetic field strengths 
at the stellar surface. Thus, according to Gil & Mitra (2001), 
such extremely strong and curved magnetic fields arc able to 
resolve in a natural way the long-standing controversy about 
the so-called 'binding energy problem', from which vacuum 
gap models suffered so far, without the necessity to involve 
bare polar cap strange stars, as proposed by Xu, Qiao & 
Zhang (1999). 

Unfortunately, only the dipolar component of the sur- 
face magnetic field, more exactly, only its component per- 
pendicular to the line of sight, is inferred from the observed 
spin-down rate of a pulsar, thus leaving the possibility for 
different hypotheses. For instance, Krolik (1991) assumes 
that the rotational energy loss of a pulsar IQCI results from 
the electromagnetic radiation generated by the ensemble of 
individual magnetic multipoles which are dominant at the 
pulsar light cylinder: this allows the surface value of each 
multipolar component to be constrained, so yielding lim- 
its on magnetic moments of high order; the tightest limits 
are obtained for millisecond pulsars, due to their short pe- 
riods P and the small values of their derivatives P. Such 
restrictions on high-order field strengths at the light cylin- 
der do not necessaxily mean their absence close to the stellar 
surface, where they can even play a dominant role. Krolik 
(1991) suggests that, if such is the case, the relative com- 



plexity of the profiles of millisecond pulsars results from the 
complexity of their actual magnetic field. 

Recent interpretations of the observations of soft X- 
rays, probably emitted due to the heating of pulsar polar 
caps, lead to similar conclusions. Page & Sarmiento (1995, 
1996) for instance, show that it is impossible to obtain the 
high level of observed pulsed fractions, or to describe asym- 
metric X-ray profiles, under the assumption of a purely dipo- 
lar field: a much better fit to the observations is obtained 
by including a suitably oriented quadrupolar field compo- 
nent. On the other hand, Cheng & Zhang (1999) in their 
model suggest that a strong and complex surface magnetic 
field significantly affects the two hard thermal components, 
usually present in the X-ray emission of pulsars. 

Interest in a multipole field dominant close to the 
star goes back to Gil (1985), who introduced a dominant 
quadrupolar magnetic field in order to explain peculiarities 
in the radio emission of the millisecond PSR 1937 + 214. 

Recently observed features of radio emission of the pul- 
sar PSR B0943-I-10, reported by Deshpande & Rankin (1999, 
2001), clearly indicate the existence of a system of 20 rotat- 
ing sub-beams, in circulation around the magnetic axis of 
the star. This well organized and stable pattern of emission 
has been observed at three different frequencies, namely at 
430 MHz, 111.5 MHz (at Arecibo) and 34 MHz (in India), 
which are supposed to represent three different altitudes 
above the stellar surface, according to the Radhakrishnan 
& Cooke (1969) model. This system of emission sub-beams 
can be associated with a set of plasma columns, with their 
feet in the acceleration region. The relativistic pair plasma 
ffow in these columns is guided by the strong dipolar mag- 
netic field up to (and beyond) the altitudes corresponding to 
the radio emission generation region. At these altitudes, in- 
stabilities occurring in a particular plasma column give rise 
to a radio emission beam, independently of similar processes 
in adjacent columns. These beams form the emission cone. 
Deshpande & Rankin (2001) point out that the columns 
themselves are not drifting or rotating: it is the action in the 
acceleration region, located near the stellar surface, which 
is doing so while it feeds and maintains the plasma column 
just above it. Several processes in the individual sub-beams 
are found to be stable: circulation time, intensities, polari- 
sation characteristics and dimensions. However, as noticed 
by Deshpande & Rankin (2001), there is no clear choice be- 
tween a continuous drift effect and particle cascade sputters. 
They speculate that the memory of the sub-beam configu- 
ration may reside in either standing waves, or in heating 
of the surface (varying along the magnetic azimuth), or in 
the seeding of the plasma columns, for instance, by sparks, 
created on a non-integer multiple of the circumference and 
thus apparently drifting. 

One way of describing successfully the observed set of 
discrete emission columns in the pulsar magnetosphere is 
to assume that close to the stellar surface the total mag- 
netic field results from the superposition of a dipolar and a 
high-order multipolar components, both fields being axially 
symmetric and star-centered. Moreover, to fulfill the require- 
ments for an efficient pair creation just above the surface of 
the star, namely the possibility to have at the stellar sur- 
face a very twisted total magnetic field, one may consider 
that the strength of the assumed multipolar magnetic field 
is higher than the dipolar field strength, or at least compar- 
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rable to it, and test that the radius of curvature of resulting 
magnetic field lines is less than the stellar radius. This con- 
cerns the domain of open field lines, as the latter includes 
the region where the emission process is supposed to arise 
(at large distances from the stellar centre, r ^ _Rs). In this 
same domain, near the stellar surface, high magnetic mul- 
tipolar components are supposed to be more efficient than 
the dipolar field in determining the properties of the vac- 
uum gap and initiating the pair production process. Let us 
remark that by 'vacuum gap' is meant the domain located 
just above the surface of the star in which a non-zero com- 
ponent of the electric field in the direction parallel to the 
background magnetic field is able to accelerate charges. 

The likelihood of having a multipolar magnetic field 
with strength higher than, or comparable to the dipole field 
strength is reasonable if one takes into account the infiu- 
ence of the general-relativistic effect of inertial frame drag- 
ging on the electromagnetic field structure near a rotating 
neutron star, as suggested by Muslimov & Tsygan (1986). 
In Section 3.2 of this paper we estimate such amplification 
factors, using the formalism developed by Muslimov & Tsy- 
gan (1986), and confirm that these factors increase with the 
multipolar order I of the multipolar components involved. 
Therefore, the actual strength of the high-order multipolar 
field components near the stellar surface may appear to be 
much higher than those estimated in the flat space-time. 
On the other hand, as mentioned above, the evolution of 
multipolar magnetic fields with I < 25 due to ohmic dif- 
fusion follows the same pattern as that of the dipolar field 
(Mitra, Konar & Bhattacharya 1999). Therefore, there are 
reasonable grounds to believe that the dipolar and multipo- 
lar magnetic field components coexist and persist with their 
relative strengths. 

The above configuration can be used to explore pulsar 
radio emission theories. For instance, the pulsar emission 
cone, when associated with a stable configuration of an en- 
semble of emission columns, could possibly be traced back to 
spark discharges in the pulsar polar gap region. This could 
be done in reference to recent works involving sparks and 
non- stationary plasma in vacuum gap models (Usov 1987; 
Ursov & Usov 1988; Asseo & Melikidze 1998; Melikidze, Gil 
& Pataraya 2000; Gil & Mitra 2001). In a different way, we 
model the whole system of observed twenty sub-beams as 
an ensemble of relativistic flnite beams, flowing along the 
curved magnetic field lines of the open magnetosphere and 
immersed in external media. The two surrounding external 
media which limit these flowing sub-beams are, on one side, 
the medium that lies on closed fleld lines, and on the other 
side, the tenuous plasma present inside the hollow cone. The 
latter surrounds the pulsar magnetic axis, and is modified 
due to the presence of high-order multipolar magnetic field 
components near the stellar surface. 

We show how, in the case of millisecond or fast pul- 
sars, the geometrical characteristics of the fiowing beams in 
the radial and azimuthal directions allow us to treat each 
beam of the ensemble as isolated and immersed in external 
media. The complex dispersion relation for perturbations of 
the system is simplified, as the beams are extremely thin and 
as the propagating waves are excited close to resonance in 
the radio domain. Such waves have simultaneously an elec- 
trostatic and electromagnetic character and thus are able to 
reach an observer. Depending on the distance of the con- 



sidered region to the surface of the star, either the 'finite 
beam', or the radiative, or the two-stream instability is the 
dominant one (Asseo 1995). Moreover, whatever the altitude 
and the concerned instability, unstable waves with charac- 
teristics in agreement with observed waves can be considered 
as candidates for pulsar radio emission. In particular, close 
to the surface of the star, above the gap region in which 
the observed beams are supposed to have their origin, the 
'finite beam' instability, dominant there, may initiate the 
radio emission process in the plasma columns. Farther away 
the classical emission processes are available. 

In the case of slower pulsars, like PSR 0943-1-10, the ge- 
ometrical characteristics of the flowing beams in the radial 
and azimuthal directions prevent considering each particular 
beam of the ensemble as isolated, just above the gap region 
and up to some distance Tj^. There, all the beams forming 
the emitting region have to be treated simultaneously. This 
can be done by introducing the 3-dimensional dependence 
of the perturbed physical quantities, leading to a simple an- 
alytical treatment of the whole system of sub-beams. In be- 
tween the gap region and the region located at and above 
r^D, simultaneous dependence of the perturbations on the 
radial and azimuthal variables will complicate the character 
of the emitted waves. In effect, the 3-dimensional character 
of perturbed physical quantities results in the coupling of 
the features of excited waves in the radial and azimuthal 
directions. Beyond the distance r^^, the sub-beams of the 
ensemble can be considered as isolated, so that the resulting 
perturbations should resemble those described above for fast 
pulsars: the character of the unstable emitted waves should 
be recovered. However, the existence of the necessary tran- 
sition between the two regions, below and above r^^, and 
thus the required continuity of the wave-solutions through 
this transition zone, suggest the possibility for an apparent 
'drift' of the observed waves. 

In the present paper we examine several magnetic con- 
figurations involving dipolar and strong multipolar magnetic 
field components, in order to determine those which pro- 
vide small curvature radii, as required by standard magne- 
tospheric models. Note that our calculations only concern 
the simplified case of the so-called aligned rotator, i.e. a 
pulsar whose dipolar magnetic momentum /x is parallel, or 
anti-parallel, to the rotation axis 17, and similarly for the 
magnetic momenta of the different multipoles. However, we 
believe that most of our qualitative results also apply to 
inclined rotators. 

As a starting point we suppose that magnetic multi- 
poles of order {l,m), created by the currents flowing in the 
stellar crust (Blandford et al. 1983; Arons 1993), are su- 
perimposed on the total dipolar field. A multipole compo- 
nent with m = we refer to as 'axially symmetric and uni- 
form', and the one with a non-zero m - as 'axially symmetric 
and non- uniform'. Let us note that uniformity, and not the 
symmetry, is the main difference between these two struc- 
tures of the magnetic field near the polar cap. Indeed, while 
m = -type fields are simultaneously symmetric and uni- 
form, m 7^ -type fields are symmetric around the magnetic 
axis but non-uniform, being fragmented as shown below. In 
Section |^, we present the general formalism to describe the 
multipolar fields. We also consider the configurations com- 
posed of the global dipolar field and a single axially symmet- 
ric and uniform multipolar field, and find out that 'neutral 
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points' and 'neutral lines', i.e. places of zero magnetic field 
strength, may exist near the neutron star surface in such a 
geometry. In Section ^, we model the magnetic structure of 
complex configurations in the vicinity of the stellar surface, 
and calculate the curvature radii of resulting field lines near 
the stellar surface. We conclude that the curvature radii of 
field lines in configurations containing only axially symmet- 
ric and uniform multipolar fields are not small enough to 
ignite the pair creation process above the pulsar polar cap, 
even if contributing multipoles are of very high order. On 
the other hand, we find that configurations involving axially 
symmetric and non-uniform multipolar magnetic fields al- 
low to achieve the small curvature radii of field lines, within 
a fraction of the modified polar cap. Furthermore, the pair 
plasma flow at large altitudes from the stellar surface ap- 
pears to be fragmented into a system of isolated filaments, 
distributed symmetrically around the pulsar magnetic axis, 
around the margin of the open dipolar fiux tube. We ap- 
ply this model to PSR B0943 -I- 10 and find resemblance 
with the observed pattern of its radio emission (Deshpande 
& Rankin 1999, 2001). In Section §, we discuss radio emis- 
sion processes in the system of discrete sub-beams of dense 
plasma, surrounded by media of reduced density, in relation 
with an analysis previously proposed by Asseo (1995). This 
is done for either millisecond and fast pulsars or for slower 
pulsars like PSR B0943 -I- 10. A summary of our results is 
presented in Section |^ 



2 MULTIPOLAR MAGNETIC FIELD 

COMPONENTS AROUND THE NEUTRON 
STAR 

Throughout this paper we use a system of spherical coor- 
dinates (r, 0, (f)) where r is the distance measured from the 
stellar centre, 6 is the polar angle (in radians) measured 
from the z-axis, and (j) is the azimuthal angle (in radians) 
measured from an arbitrary origin. The z-axis is directed 
along the dipolar momentum of the star /x. 



2.1 General description of multipolar fields at 
close distances from the stellar surface 

In the near zone limit {kr <^ 1), the magnetic field vec- 
tor, B^"^{r,6, 4>), associated with a given magnetic multi- 
pole (Im), can be expressed in terms of spherical harmonics 
Yim{e,4>) (Jackson 1975): 



B''"(r,e,< 



„! + l 



(1) 



In a spherical geometry, the components of an in- 
dividual multipolar magnetic field vector B'"^ (r,9, (j>) = 
{B';"ir,e,^),B'frir,e,cb),B';"{r,9,^)) are written: 



Bl-ir,e,cP)^-A.l±±^Y,40,cb 

Bi-{r,e,^) = ^^Y^^/\e,<p), 

47r gi„ 



B';^ir,e,4>) = 



-im sin 6Yim{6, 



(4) 
(5) 
(6) 



21 + 1 r'+2 

Particular multipolar magnetic components are ob- 
tained for specific values of I and m, with |m| < I. Namely, 
dipolar components are obtained for / = 1 and either m — 
or m = ±1, quadrupolar components for I = 2 and either 
m — 0, m = ±1 or m = ±2 and so on. If m = and I is 
arbitrary, then one obtains axially symmetric and uniform 
multipolar components. 

The magnetic moment, associated with the strength of 
the magnetic multipole component (Im) at the surface of 
the star, is defined as = Bi'''B}+'^. 

The total magnetic field of the neutron star can be writ- 
ten as a sum of the individual multipolar magnetic field 
components. In the hypothesis where all the multipoles are 
centered at the centre of the star, 



B{r,e, 



E 



E 



Sim. / 



(7) 



i = l m=0,±l,...±i 



Magnetic field lines for the total magnetic field are ob- 
tained, as usual, from an integration of the set of differential 
equations: 

dr rd9 r sin 9d(j> 



Br{r,e, 



Be{r,e,cf,) B4r,9,q>)- 



(8) 



A unit vector in the direction of the local magnetic field 
B(r,9,(j)) can be defined as: 



b{T,9,^) = 



B(r,9, 



(9) 



B{r,9,cj>) 

The absolute value of the derivative of b(r, 9, (f>) along 
the field line is nothing else than the curvature of the mag- 
netic field line. At a given point of space (r, 9, (f)), the radius 
of curvature of the total magnetic field line is therefore de- 
fined as 

1 



p{r,9,cb) 



\{b{r,9,q>)V)b{r,9,<j,)\- 



(10) 



Naturally, if one assumes the existence of only one pure 
individual multipolar magnetic field component around the 
star, the equations of its field lines and their curvature radii 
can be easily obtained from equations (^ and (|l^) by sub- 
stituting B''"(r, 9, (p) instead of the total field B{r, 9, (j)). 



Spherical harmonics are written in terms of the associated 
Legendre functions, as 



(9, 0) = ^ |_^Pr (COS 9) e-^ 



(2) 



where (x), the associated Legendre function, is defined 
as 

V2 d'+'" 



d2;^^ 



fx^-lV. 



(3) 



2.2 The case of axially symmetric and uniform 
multipolar magnetic fields: magnetic field 
components and the magnetic vector potential 

In this section we treat multipolar magnetic fields of type 
Bio, i-e. those for which I is arbitrary but m = Q. Such 
fields can be regarded to as 'dipolar-like', as they are axially 
symmetric and uniform in the vicinity of pulsar polar caps, 
like the dipolar magnetic field. 
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According to equations — Q), for the case with m = 
0, the components of the total magnetic field close to the 
surface of the star can be expressed in a general form as a 
sum of multipolar magnetic field components, namely 



oo 



oo /2n<l 

-J^ f [2^ C,_2„ COS^ 



(11) 



¥°{r,e) = rsmeA'^{r,9), 



(16) 



which is an integral of the field equations. Effectively, using 
the fact that magnetic field components of the multipolar 
magnetic field with indices (10) derive from a vector poten- 
tial A''^{r,9), we obtain: 



1 



r sm a 

_ld_ 
r dr 



^K°(r,e)sine) 



^y{r,e),(n) 



1 



- sin d: 



^*"'(r,e), (18) 



Be{r,e) = Y,Bf{r,e) = 



E 



c /i sm f) iQ 



^2n + l<! 



^1 E 



Ci-l-2n cos 



{i-l-2n) 



B4r,e) = ^B;°(r,e) = 0, 



(12) 



(13) 



where /x'" — BI'^RI'^^ is the magnetic moment associated 
with the strength of the multipolar magnetic field com- 
ponent at the surface of the star, /J^ being a positive quan- 
tity according to its definition; is a constant which de- 
pends on the order I of the multipole; and F'^ are specific 
functions which only depend on powers of cos 9; c;_2n and 
c;_i_2n are numerical coefficients. 

As implied by Maxwell's equations, the static magnetic 
field vector can be defined through the magnetic vector po- 
tential A, 8S B — {V X A) . Thus, one can associate a vector 
potential A'"^ with any arbitrary magnetic multipole with 
indices I and m. Obviously, in the axially symmetric and 
uniform case with arbitrary I and with m = 0, each mul- 
tipolar component of the magnetic field vector, B'''\r,9), 
derives from a potential vector A'°{r, 0). Assuming that the 
magnetic field is created by an extremely thin current that 
circulates all over the crust of the neutron star, the potential 
vector has only an azimuthal component A^{r, 9), which can 
be specified for each particular magnetic multipole having 
indices I > 1 and m = 0, in terms of the corresponding cur- 
rent, or equivalently in terms of the magnetic moment 
and the variables 9 and r, as 



aU), o\ M sm9 10 



^2n + l<l 

E 



Cl_i_2n cos 



(i-l-2n) 



.(14) 



Bf{r,9)^0. 



(19) 



The function Vl/"'(r, 6^) is constant along a given field 
line; the behaviour of multipolar magnetic field components 
B';9{r,9) and B^g{r,9), necessarily follows the behaviour of 
the derivatives 9r*'°(r, 9) and de'^'°ir, 9) along a given field 
line and in a given plane in which (j> = const. As a conse- 
quence. 



¥°(r,9) = const 



(20) 



gives the general form for multipolar magnetic field line 
equations. For a magnetic multipole (/O), from equa- 
tions (jlj) and (p^, equation ( |2c| ) can be written 



,n o /2n + l<i 

— — ^ I 



^ C;_i_2„COS^ '\ 



— const. (21) 



The functions F are very easily obtained for the suc- 



cessive multipolar components: F 



1, F^ 



cos 9, F-" 



3 + 5 cos 2^, . . . , F'° = E^lV^' Ci-i-2n cos (/ - 1 - 2n)9, 
making use of the transformations 

of circular functions F"^(^J2l"Jo'^' c,_i_2„ cos^'-^-^") 9^ = 

P'''°(j2T=o^^' '=;-i-2n cos {l-l- 2n)9^ . These functions, 
according to equation (|2^), provide magnetic field line equa- 
tions for the appropriate particular multipolar fields. Draw- 
ing of the field lines in a plane where <j) ~ const, shows that 
a pure magnetic multipole of order / makes 21 loops around 
the centre of the star. 

A function 5'(r, 9) can be associated with the total mag- 
netic field: considering that the total magnetic field, in a 
given plane = const, is the sum of the successive multi- 
polar magnetic field components, the function ^'(r, 0) is the 
sum of the successive functions ^'"'(r, 6') defined above. The 
general form of the function ^'(r, 6') can be written 



where the functions F"^ have been introduced just above. 

Magnetic field lines are obtained as usual from an in- 
tegration of the set of two differential equations (^ . Here, 
for the ajcially symmetric and uniform case with arbitrary 
I and m = 0, the (^-component of the magnetic field being 
zero, it is sufficient to consider the multipole magnetic field 
components in a plane (j> = const, and to integrate only the 
first differential equation, namely 



dr 



rd9 



B'°{r,9) B'°{r,9)' 



(15) 



A systematic derivation of the equations of field lines 
of any particular magnetic multipole can also be obtained, 
introducing the function 



*(r,e») = ^*"'(r,e) =rsm9"^A'^\r, 

1=1 1=1 



(22) 



Consequently, the total function '^{r,9), written as an infi- 
nite sum. 



*(r,0) = sin^e^^F'° 




C;_i_2„ COS 



(i-l-2n) 



,(23) 



provides the total magnetic field line equation '^{r,9) = 
const in the given plane (p = const. 

Using the specific contributions of the first multipoles, 
together with the functions introduced in equation ( |l2| ) 
(see also the text below equation Ell) , and with the adequate 
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numerical coefficients, we have: 



^ + iij- COS 6» + ^^3- (3 + 5 COS 26*) + 



*(r, 6*) = sin 

+ ^^"1 c;_i-2„cos(/-l-2n)e) + 



/2n + l<i 



2.3 The case of axially symmetric and uniform 
multipolar magnetic fields: neutral lines, 
neutral points 

The surface 'ii{r,9) — separates field lines with opposite 
direction, that is field lines on which the magnetic field has 
reversals. Whatever the number of multipoles, the function 
'i'{r,9) is zero for either 6 = 0, or 6 — n, or when the 

variables r and 6 axe such that J2Zi'^io'"^(^y^) = 0- Ob- 
viously, there should be some particular relation between 
the magnetic moments /i'" to have 'i>{r,6) = 0, and thus a 
total magnetic field with zero strength, B{r,6) — 0, in any 
arbitrary plane = const and at some point {r,9), as, for 
instance, at the surface of the star, where r = R^. 

Values of the magnetic field on the lines 9 — 0, or 9 = n 
are easily obtained. From the equations (jl^) - (jl^), giving 
the magnetic field components, it is clear that on the lines 
9 = and 9 = -n, B^{r,9) = 0, as it is zero everywhere, 
Be{r,9) = 0, since sm9 is in factor, but Br(r,9) takes-on 
values that depend on the particular point where the ra- 
dial component of the multipolar magnetic field is calcu- 
lated. In particular, on the line 9 = n, there is a point at 
which Br{r,9) = 0. A study of the behaviour of the mag- 
netic field around this point will probably show that it may 
be considered as a neutral point for particular configura- 
tions of the magnetic field in the pulsar magnetosphere. For 
instance, this should be the case when the total magnetic 
field can be represented by a combination of either dipolar 
and quadrupolar magnetic fields, or dipolar and octupolar 
magnetic field, or dipolar and 16-polar magnetic field and 
so on, the equations for the field lines being, respectively. 



^{r,9) = sin 



■ cos 9 



+ ^(3 + 5 cos 261) 



*(r,6') = sin 9 I — + ^ (9 cos -f 7 cos 36') 



(25) 



(26) 



(27) 



and so on. 

On the side of the pole where 9 = 0, there is no 
real solution for r, the r-solutions being either negative 
or imaginary. By contrast, on the side of the other pole, 
where 9 = tt, there is a point r at which the magnetic 
field strength is zero, and changes sign around this point. 
Effectively, for the combination of a dipole plus a multi- 
pole of even order (/ = 2n,n = 1,2,...), there is a point 
where the strength of the magnetic field is zero. The solution 



lip = (numerical coefficient) (^Bl° /B. 



10\l/(!-l) 



Rs, depends 



on the ratio of the multipolar and the dipolar magnetic field 
strengths. More generally, for an arbitrary 6-angle there is a 
neutral line on which the magnetic field is zero. For instance. 



if the magnetic field is taken as the simple combination of 
dipolar and quadrupolar magnetic fields, the function 



*(r,e) 



+ 



(28) 



(24) becomes zero at either ^ = or 6 = tt, or on the line where 



+ 



= 0, (29) 
that is to say, on the circle defined by the equation 



-El. 

^2 



r 

r'.W' 



(30) 



In this case, there are three neutral lines: the two axes, 
9 = and 9 = n and a circle in the plane {r,9), on which 
the only acceptable r-solutions correspond to S-angles in 
the domain [7r/2, 37r/2]. Indeed, depending on the respective 
values of fiio and /i20, the r-solutions on the neutral circle 
have a physical meaning, or none. As an example, for the 
combination of dipolar and quadrupolar magnetic fields, the 
'neutral point' is located at distances from the centre of 
the star rll = 1.16 (B^^/BI") R^. On this line 9 = n, Br 
is positive at smaller distances, r < r^p, and negative at 
larger distances, r > r^p. The 'neutral point' is located at 
the surface of the star only if ° is less than " , so that 
Bf = _B,^°/1.16 = OmB^°. 

As clear from the general results given above, similar 
locations for the 'neutral points' can be easily obtained for 
any combination of a dipolar field plus a multipolar field 
with an even index A study of the behaviour of the mag- 
netic field on lines close to the 'neutral line' is necessary to 
conclude about the specific physical processes in this region. 



3 MODELLING THE MAGNETIC FIELD 
STRUCTURE NEAR THE STELLAR 
SURFACE: TOPOLOGY OF POSSIBLE 
CONFIGURATIONS AND CURVATURE 
RADII OF FIELD LINES 

3.1 General description of the method of 
calculation 

In this section we study the topology of magnetic field con- 
figurations near the stellar surface, which consist of a multi- 
polar magnetic field ij''" (with I > 1 and |m| < /), superim- 
posed on the global dipolar field of the pulsar B'^ = B^°. We 
also introduce the ratio of the total and the dipolar magnetic 
field strengths. 



/3(r,e,0) 



B{r,9,^) 



Sd(r,0, 



(31) 



where B{r,9, 



\B'^{r,9,(t)) + B''"(r, 0,0)1. Therefore, 



f5{r,9,(f)) indicates how many times the strength of the to- 
tal field exceeds the strength of its dipolar component, at a 
given po int of sp ace (r, 9, (f)). In the calculations reported in 
Sections 



3.3 and 3.4 



particular attention is paid to the cases 
where the curvature radius of magnetic field lines at the stel- 
lar surface is smaller than the stellar radius, i.e. ps < Rs, 
at least within a fraction of the modified polar cap. Intu- 
itively, one would expect that such configurations assume 
a significant contribution of the multipolar component, so 
thaX I3s = P{Rs, 9, 4)) > 1. 
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In order to test various magnetic configurations, we re- 
define tlie expression for tlie multipolar field components 
(equation ^ by introducing an extra amplification factor ^ 
as a free parameter, so that this expression now takes the 
following form: 



(32) 



A specific choice of the ^-factor in every particular example 
considered below, is justified by the output: it should al- 
low us to obtain the expected structure and strength of the 
surface magnetic field above the modified polar cap. Am- 
plification of certain multipolar components is discussed in 
detail in Section ^.2| . 

In order to model the geometry of the total magnetic 
field, one has to solve the set of differential equations 
which is now convenient to rewrite in the following form; 



de(r) ^ Be (r, e(r), 0(r)) 
dr Br{r,e{r),cl,{r))' 

d0(r) _ B4r,e(r),0(r)) 



r sin S(r) 



dr 



B,.{r,9(r),cbir))' 



(33) 
(34) 



This set of equations has to be solved for the unknown func- 
tions 9{r) and (;/>(r). 

Of course, the complex surface magnetic field should 
transform into the global pure dipolar field at large alti- 
tudes from the stellar surface, where multipolar components 
effectively vanish. This is warranted by the choice of correct 
boundary conditions while solving the set of equations ( |33| ) 
— (^) . Since the generation of radio emission is associated 
with the bundle of the open dipolar magnetic field lines (i.e. 
those crossing the light cylinder), we restrict our calculations 
to this domain of the magnetosphere. 

Let us introduce some distance r^, where sur- 
face fields are negligibly smaU, B''"(r,, 6i(n), (?!)(n)) < 
B'^{r^,9{r\),(l>{r\)), so that the total magnetic field can be 
treated as purely dipolar there. At this distance, the points 
of the last dipolar open field lines delineate a certain circle, 
and their azimuthal angle (/"(rj takes arbitrary values within 
the range < (/"(rj < 2tt. According to the equation for the 
dipolar field lines. 



const, 



the polar angle of these points can be written 



(35) 



arcsin 1.45 x 10" P" 



Here R^^^ — Pc/2n is the light cylinder radius of the pulsar, 
and P is the pulsar period. The dipolar field falls off with 
the distance as (Ra/r)^ , whereas the multipolar field of order 
I falls off as (Ps/r)'"*"^, according to equations (^ — (j^. 
Therefore, using equation (^l|), the characteristic distance 
r^, where S'™ ~ B'^, can roughly be estimated as 



(/3s 



Rs 



(37) 



Thus, the dipolar field dominates over the multipolar field in 
the domain of altitudes r ^ (and vice versa). Choice of r. 
in the latter domain warrants a smooth transition from the 
complex surface field to the open dipolar field lines. Hence, 
this is equivalent to the integration of the field line equations 
from the stellar surface up to the light cylinder. 



As well-known, observations of pulsar radio emission 
and most of the models for its generation imply that it orig- 
inates at altitudes r > 50i?s. Interpretation of radio obser- 
vations indicates that the magnetic field is purely dipolar in 
the generation region, so that complex surface fields should 
vanish at such altitudes^. In other words, the complex sur- 
face field lines should connect to pure dipolar field lines well 
below r = 50-Rs. We verify that this requirement is fulfilled 
through out our calculations, for each chosen ^-parameter. 

Hence, the boundary conditions at r = Tj for different 
field lines are to be taken within the following angular range: 



< 0(r.) < e., 
< 4>(r.) < 2n. 



(38) 



Each set of boundary conditions defines a single open field 
line, starting at the stellar surface, where the multipolar 
magnetic field dominates over the dipolar field, and extend- 
ing to the domain located at altitudes r > where the field 
is purely dipolar and axially symmetric and uniform. We 
use a Fehlberg fourth-fifth order Runge-Kutta method (see, 
e. g., Forsythe et al. 1977) in order to obtain the numerical 
solutions of the set of equations ( ^3| ) - (^). 

The results of our simulations are presented in the fol- 
lowing subsections. We separate them into two cases, which 
have diverse physical consequences, namely: i) configura- 
tions composed of axially symmetric and uniform multipolar 
magnetic fields with m — and the global dipolar field, and 
ii) configurations composed of axially symmetric and non- 
uniform multipolar fields with m ^ and the global dipolar 
field. Before starting to discuss these particular cases, let us 
present one of the arguments in favour of an anomalous am- 
plification of the 'weight' of multipolar components at the 
stellar surface, with respect to the dipole field. 



3.2 Intensification of multipolar magnetic field 
components due to the general relativistic 
effect of inertial frame dragging 

As mentioned in Section ^ high-order multipolar fields may 
be intensified close to the surface of a rotating neutron star, 
relatively to their value in a fiat space, due to the general- 
relativistic effect of inertial frame dragging studied by Mus- 
limov & Tsygan (1986). The exterior gravitational field of 
a neutron star, rotating with a constant angular velocity fi, 
is described by the metric, which is locally inertial at the 
infinity: 



ds^ = gai3dx°'dx'^ , a, /3 = 0, 1, 2, 3. 



(39) 



The nonzero components of the metric tensor are written as 
follows: 



500 



-gii = 1 



Rei 



922 



2 

-r , 



2 ■ 2 a 2GJ . 2 o 

533 = -r sm 9, gos = sm 9. 



(40) 



* Although, the anomalous dispersion measure of the high- 
frequency versus low-frequency radio pulses, observed in some 
pulsars (Davies et al. 1984; Kuzmin 1992), has been attributed to 
the presence of the quadrupolar component of the magnetic field 
in the emission region. 
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Figure 1. Amplification of different multipolar components due 
to the frame dragging effect (equations H and||). The values of 
the functions fi{r) and ,/annai(r) at the stellar surface (r = Rb), 
given by equations ( ^ ) and (^) , are plotted versus I in panels (a) 
and (b), respectively. This plot corresponds to the neutron star 
with the presumed mass M = IAMq. The curves with squares, 
circles and triangles represent the calculations for the presumed 
stellar radii of 7 X 10^ cm, 10^ cm and 1.3 X 10^ cm, respectively. 



Here Rg = 2GM/c? is the Schwarzschild radius, G is the 
gravitational constant, c is the speed of light, M is the mass 
of the neutron star and J is its angular momentum. 

Muslimov & Tsygan (1986) solved the Maxwell equa- 
tions in such a geometry, and found the modified magnetic 
field components modified as follows: 

B,.(r,e,0) = y^(Z + l) (r^y'^ViW himYim{e,\), (41) 



Im 



Im 



(^) 



1+2 



/googi{r) 6;™ Ve,0Yi™(6', A). (42) 



Here are some coefficients defined from the boundary 
conditions (see equations |^ - A = <?!> — 0,t, and the func- 
tions fi{r) and gi{r) can be expressed in terms of hyper- 
geometric functions, 

Mr)=F(^l,l + 2-2{l + iy,^Y (43) 



(44) 



gi{r)=F(^l + l,l + 2;2{l + l)-^y 

Therefore, the modified radial and angular multipolar 
magnetic field components are connected with the corre- 
sponding components in the flat space, as 



Br ^ fi{r)Br , 
e,4> = V9oo9i(r)Bg^^. 



(45) 
(46) 



Estimates of the values of the functions fi{r) and 
/Qmgiij') at the stellar surface, r = Rs, are displayed in 
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Figure 2. The same as in the previous figure, but for the neutron 
star with the presumed mass M = 1.6Mq. 



Figs. |l| and |^ for a presumed mass of the neutron star 
M = 1.4Mq and M = 1.6Mq, respectively. The functions 
ft{r) and gi{r) are calculated for three different cases of 
the model neutron star radii, namely Rs = 7 x 10^ cm, 
10^ cm and 1.3 x 10® cm. First, one notices that the re- 
sulting intensification is about the same for the radial and 
angular magnetic field components, and is stronger for the 
higher-order multipoles. Second, multipolar fields are ampli- 
fied more intensively near the surface of relatively smaller 
and more massive neutron stars. Indeed, according to Fig. ^ 
high-order multipolar magnetic field components with I < 25 
in the curved space may be even 10* — 10^ stronger than the 
corresponding components in the fiat space. 

The evolution with time of multipolar magnetic field 
strengths in isolated neutron stars has been studied by Mitra 
et al. (1999), assuming that the field is mostly generated in 
the outer crust of the star and its later evolution is due 
to the ohmic decay of currents in the crustal layers. They 
have shown that if one assumes the same initial strength 
for all multipolar magnetic fields, the evolution is such that 
the reduction with time is similar to that of the dipole field 
component, except for very high multipole orders (/ > 25). 

A combination of the two described effects, that is, the 
intensification of multipolar magnetic field components due 
to general relativistic effects and evolution with time of their 
strengths, suggests that multipolar magnetic field compo- 
nents with high order should be present at the surface of 
neutron stars and important enough to determine the mag- 
netic configuration in their vicinity. This justifies the choice 
of ^ (equation ^) as a free parameter in our calculations 
below. 
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3.3 Complex configurations containing axially 

symmetric and uniform multipolar magnetic 
fields: field line topology and curvature radii 

In this section we present the results of our calculations for 
the magnetic configurations composed of the global dipolar 
field and a surface multipolar field with m = Q and arb itrary 
Z > 1. As mentioned in the beginning of Section 2.2, such 
fields can be regarded as 'dipolar-like', as they are axially 
symmetric and uniform in the vicinity of pulsar polar caps. 
Due to this fact, the ratio of the total and the dipolar r nag - 
netic field strengths at the stellar surface, (3s (equation 
is independent of the azimuth and also almost indepen- 
dent of the polar angle 6, within the modified polar cap of 
a pulsar. It should be noticed that the polar cap still has 
a circular shape, as in the pure dipolar case. However, as a 
result of magnetic flux conservation, the presence of a strong 
surface magnetic field with /3s > 1 reduces its angular radius 
to = 6p[3~'^'^ , compared to that of the pure dipolar polar 
cap, 6*^ « 1.45 X lQ-'^p-°-^ [rad]. 

For pure multipolar magnetic field lines with m = 
(i.e., in the absence of the dipolar component), a first-order 
approximation for the curvature radius is obtained from 
equation (|lO|), as 



(47) 



/(/ + 2) 



For example, the curvature radius of the last open field lines 
at the stellar surface {r = Rs,9 = 9^p) for a pure multipole 
{10) is: 



1(1 + 2) 



-Rs 



(48) 



However, as demonstrated below, the presence of the dipolar 
magnetic field in the actual configuration near the stellar 
surface modifies the values given by equation (^), especially 
assuming small and moderate values of /3s- 

As an example, we consider a sample pulsar with P = 
0.2 s and P = 10~^^ s/s. The surface dipolar magnetic field 
of such a pulsar is estimated as ~ 3.2 x 10^® i^^) ~ 
4.4 X lO^'^ G. Such a value allows intensification of the sur- 
face field up to 100 times, without exceeding the critical field 
Bq = 4.4 X 10^^ G. We model the dipole plus single mul- 
tipole magnetic field lines of t his pulsar, according to the 
scheme described in Section 



3.1 



and calculate the curvature 
radii of the last open field lines at the stellar surface. Our 
aim is to check how strong should be the contribution of 
axially symmetric and uniform multipolar components near 
the stellar surface, in order to achieve the desired small val- 
ues of curvature radii of the resulting open magnetic field 
lines, as required for an efficient pair creation process. 

The curvature radii of the last open field lines at the 
stellar surface are presented in Figs, ^(a), (b) and (c), which 
correspond to cases where the total surface magnetic field 
exceeds the dipolar component 2, 10 and 100 times, respec- 
tively. The solid lines with circles represent the results of 
our numerical calculations for given complex configurations. 
The dotted lines with triangles show, for comparison, what 
would be the curvature radii of the involved multipolar mag- 
netic field lines in the absence of the dipolar field (as given 
by equation ksl) . 

In Fig. ffa) it is seen that, if the contribution of surface 
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Figure 3. Curvature radii of the last open field lines at the stellar 
surface in the configurations composed of the global dipolar field 
and the axially symmetric and uniform multipolar field with m = 
0, as functions of the multipolar order I. The three panels (a), 
(b) and (c) correspond to the three dilferent cases of relative 
contribution of the multipolar components to the surface field 
strength, /3s = 2, f3s = 10 and /3s = 100, respectively. The solid 
lines with circles correspond to the calculations for the complex 
surface field, whereas the dotted lines with triangles represent 
curvature radii, as they would be in the absence of the dipolar 
field (equation |48[). 



multipolar fields is comparable to the dipolar one, /3s = 2, 
the radius of curvature of resulting field lines remains larger 
than the stellar radius, pa > Rs, even for multipolar fields of 
very high order, I > 25. On the other hand, the dotted line 
shows that, in the absence of the dipolar field, ps < Rs would 
be achieved already for multipoles with I > 12. Therefore in 
this case, multipolar fields are simply too weak to 'twist' the 
field lines strongly enough. 

The situation is somewhat better, if one assumes a 
stronger contribution of multipolar fields, /3s — 10. This case 
is presented in Fig. ^b). As seen in this figure, the curva- 
ture radii of field lines become less than the stellar radius, if 
the multipolar order of contributing multipoles is very high, 
I > 24, still quite a possible case, according to Mitra et al. 
(1999). 

However, a further increase of the contribution of sur- 
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face multipolar fields worsens the situation. For instance, 
in the case of /3s = 100, presented in Fig. ^(c), curvature 
radii never drop below the stellar radius, even when very 
high-order multipolar components {I > 25) contribute into 
the total magnetic field. This fact is easy to explain: such 
strong multipoles fully determine the surface field structure 
near the stellar surface. Indeed, the dotted line in Fig. ^(c), 
corresponding to the absence of the dipolar component, al- 
most follows the solid line, corresponding to the total sur- 
face magnetic field. On the other hand, the high surface field 
strength leads to a highly (10 times) reduced angular radius 
of the polar cap, which, according to equation (US) , provides 
larger curvature radii. 

Therefore, we conclude that in the configurations com- 
posed of the global dipolar field and the axially symmet- 
ric and uniform multipolar fields, the curvature radii of the 
open field lines do not generally have curvature radii, small 
enough to fulfill the conditions for a copious pair creation 
process. At the same time, as seen in Fig. ^(c), in the con- 
figurations composed of global dipole and axially symmetric 
and uniform multipole, it is impossible to have both ps < Rs 
and /3s ~ 100 in the same configuration, as implied by the 
vacuum gap model of Gil & Mitra (2001). These conclusions 
are true at least for multipolar fields of order I < 25, which 
do not significantly dissipate during the pulsar characteristic 
life-time (Mitra et al. 1999). 



3.4 Complex configurations containing axially 
symmetric and non-uniform multipolar 
magnetic fields: field line topology and 
curvature radii 

Let us now study magnetic field configurations composed of 
a global dipolar field on which non-axisymmetric multipolar 
magnetic fields (i.e. those with m 7^ 0) are superimposed. 
We again restrict our calculations to the open magnetic fiux 
tube. As in the previous section, we start our calculations 



at some altitude r. where B ' 



« B^(n 



and 



solve the set of equations (|3^) and ( |34| ) from down to the 
stellar surface. While doing so, we take the boundary condi- 
tions in the range (^) . This allows us to find the footprints 
of the tube of open field lines on the stellar surface, i.e., to 
construct the modified polar cap. Our aim is to get solutions 
which provide p < _Rs at the stellar surface, at least within a 
fraction of the modified polar cap, which is necessary for an 
efficient pair production process. We adjust the free param- 
eter ^ (introduced in equation ^) in each particular case 
considered below, in order to obtain desired configurations. 
Throughout our calculations we control the value of the total 
surface magnetic field, so that it does not exceed Bq within 
the polar cap. Obviously, there is a great variety of possible 
configurations. Here we present the most illustrative ones. 
Let us first make a few general remarks. 

A pure multipolar magnetic field {l,m) with m 7^ 
has 2m(Z — m poles over the whole stellar surface. Out 
of those, along each latitude containing poles, 2m poles of 
alternating polarity are distributed around the axis of sym- 
metry. At a given latitude, each pole with a given polarity 
is neighboured on its sides by two poles of the opposite po- 
larity. The inclusion of the global dipolar field changes this 
topology. In effect, the dipolar field lines which enter one of 
the two polar regions on the stellar surface, are only able 




Q.2 

01 the peripheral part of the 



Figure 4. Magnetic field structure 1 
open flux tube near the stellar surface, corresponding to the mix- 
ture of the dipolar field and the multipole of order Z = 4 and 
m = 4. The Cartesian coordinates x, y and 2 are represented in 
the units of the stellar radius R^. 
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Figure 5. Peripheral fraction of the modified polar cap for the 
contributing multipole Z = 4, m = 4. The footprints of the open 
field lines are represented by circles. Deformation of the polar cap 
under the infiuence of the local multipolar poles is clearly seen. 



to connect to the local multipolar poles of a given polarity, 
and not to the poles of the opposite polarity. This reduces to 
m the number of local poles, which are close to each of the 
polar regions. One expects that the influence of these local 
multipolar poles is stronger at the annular edge of the polar 
cap. Consequently, the polar cap becomes modifled: its edge 
becomes rather fragmented into m independent 'hot spots', 
as demonstrated below. 

The most interesting cases, in our opinion, are those 
where the maximum number of local poles surround the 
magnetic z-axis of symmetry of the dipolar fleld at the clos- 
est angular distance, i.e. those with Z = m. As an example, 
let us first consider the multipolar magnetic field with Z = 4 
and m = 4, superimposed on the global dipolar field. 

The open field line geometry for such a configuration 
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near the stellar surface and the corresponding modified po- 
lar cap (composed of the footprints of the open field lines) 
are presented in Figs. ^ andj^, respectively. The calcula- 
tions were done for a sample pulsar with P = 1 s and 
P = 10~^^ s/s. For simplicity, we only present in these fig- 
ures the outermost 0.2 fraction of the open flux tube and 
its projection at the stellar surface (the inner, more cen- 
tral open field lines are less affected by the multipolar field, 
and follow rather the dipolar pattern). In other words, the 
boundary conditions at an altitude = lOiis (see equa- 
tions ^ and liil ) were taken within the following angular 
ranges: 



< 0(r-i) < 27r, QM, < e{r.) < 



(49) 



Each of these ranges was divided into a fixed number of 
equal parts, both by the azimuth and the polar angle, so 
that the angular spacings between the initial points at the 
altitude r. were taken equal. 

Fig. ^ shows that the bundle of last open field lines, 
which is circumferentially uniform at higher altitudes, r > 
2Rs, where the multipolar magnetic field is negligible, tends 
to be deformed under the influence of the multipolar field 
closer to the stellar surface. As it is also seen in Fig. ^ the 
field lines are most significantly deviated outwards in those 
0- direct ions, where the poles with an appropriate polarity 
are located. On the contrary, in the directions towards the 
multipolar poles of the opposite polarity, field lines deviate 
inwards from the pure dipolar circular pattern. The resulting 
polar cap appears to be rather distorted, as displayed in 
Fig. I 

One would expect that both the magnetic field strength 
and the curvature radius of the field lines also vary within 
the periphery of the modified polar cap, so that both of them 
now appear to be functions of and 0, as opposed to the con- 
figurations containing axially symmetric and uniform multi- 
polar components (discussed in Section 3.3), where both I3s 
and ps could be treated as constant throughout the modified 
polar cap. Let us study this dependence. 

Fig. ^a) displays the azimuthal variation of /3b (i.e., the 
value of the total surface magnetic field normalized to the 
local value of the dipolar component) at the feet of the last 
open field lines. Fig. ^(b) displays the same total magnetic 
field measured in the units of the critical field _Bq. An appar- 
ent increase of the magnetic field is seen in the direction of 
the appropriate multipolar poles, whereas the field retains 
its almost dipolar strength in between these peaks (the min- 
imum field strength is achieved towards the multipolar poles 
of the opposite polarity). In Fig. ^(b) we see that the field 
may increase to a significant fraction of Bq on certain field 
lines. 

Fig. ^(c) shows the azimuthal variation of the curva- 
ture radius (derived from equation [lo[ ) of the last open field 
lines at the stellar surface. Clearly, ps ~ (0.1 1)-Rs only 
for the group of field lines stretched in the direction of the 
multipolar poles, whereas the surface value of the curvature 
radius of the remaining field lines (directed rather to oppo- 
site poles) is much larger than the stellar radius. The same 
conclusion would follow from an intuitive analysis of Fig. ^ 
where it is seen that the multipolar poles of an 'appropriate' 
polarity tend to bend the field lines in their direction, thus 
to increase their curvature. On the contrary, the opposite 
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Figure 6. Dependence of Ps, Bs/Bq and ps on the azimuthal 
angle (j>, at the feet of the last open field lines, for the contributing 
multipole Z = 4, m = 4. Conditions for pair creation arc only 
fulfilled in localised regions, where ps < Rs- 



poles 'repel' and thus straighten the field lines directed to 
them. 

Dependence of /3s, B^/Bq and ps on the polar angle 9 
in the direction of an 'appropriate' multipolar pole is dis- 
played in Figs. I^a), (b) and (c), respectively. These figures 
show that l3s > 1 and ps ~ (0.1 -t- 1)Rs at the peripheral 
fraction of the polar cap. At the same time, as we have seen 
in Figs, ^a), (b) and (c), the polar cap periphery is also 
modulated azimuthally. 

Therefore, it can be concluded that the conditions for 
a copious pair creation process are only fulfilled in local- 
ized regions, symmetrically distributed around the dipolar 
magnetic axis, within the thin peripheral part of the modi- 
fied polar cap. Although at high altitudes above the stellar 
surface the field lines strictly follow the dipolar geometry, 
which is azimuthally uniform, such a non-uniform pattern 
of plasma creation near the stellar surface is to be reflected 
in the non-uniformity of the secondary plasma density dis- 
tribution within the cross-section of the open magnetic flux 
tube at some large altitude (r > r.). First of all, the pair 
plasma will only flow along the outermost narrow conical 
layer of the open dipolar flux tube. The radial S-width of 
this layer can be estimated from Fig. ^(c). Indeed, out of 
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Figure 7. Dependence of /3s, Bn/Bq and pa on the polar angle 9, 
in the direction of the local multipolar pole of order I = A,m = A. 



the 30 footprints of the open field lines displayed in this 
figure, only outermost 10 of them are feet of the field lines 
with ps ~ (0.1-hl)i?s. Recalling that these points have equal 
spacing at large altitudes and also taking into account the 
boundary conditions (^^, one can estimate that the width 
of the conical layer loaded by the pair plasma will consti- 
tute about (1 — 0.8)/3 ~ 0.07 fraction of the total angular 
radius of the open flux tube. Moreover, this flow will also 
be fragmented azimuthally: pair plasma will flow only along 
the bundles of field lines which were stretched towards the 
local multipolar poles while emerging from the stellar sur- 
face. In the case of ? = 4, m = 4, presented above, this 
corresponds to four isolated bundles of field lines. There- 
fore, plasma flow will be confined into a thin 'patchy hollow 
cone'. On the other hand, it is well known that physical pro- 
cesses leading to the generation of pulsar radio emission are 
associated with the existence of a dense pair plasma. That 
being so, the resulting emission cone will appear as split 
into symmetrically distributed 'hot spots'. This is how, in 
our opinion, the 'memory' of the surface field structure is 
preserved in the pulsar radio emission beam. 

Let us notice that such a geometry is able to provide 
azimuthally localized emission beams in both the vacuum 
gap model (Ruderman & Sutherland 1975) and its com- 
peting space charge limited flow model (Arons & Scharle- 



mann 1979, modified by Muslimov & Tsygan 1992). As 
well-known, the latter model has previously been criticized 
mainly because of its inability to provide isolated plasma 
filaments and associated beams of emission. The inclusion 
of multipolar magnetic fields with m 7^ removes the ne- 
cessity of introducing spark discharges above the polar cap, 
as complex surface fields directly fragment the fiow into iso- 
lated plasma columns. 



3.5 Complex configurations containing axially 
symmetric and non-uniform multipolar 
magnetic fields: a model for PSR B0943 + 10 

A study of field line topology and determination of curvature 
radii in a configuration involving a multipole of high-order 
is proposed with intent to interpret the special radio emis- 
sion features observed for PSR B0943-I-10, whose dynamical 
parameters are: P = 1.0977 s and P = 3.529 x 10"^^ s/s. 
We assume that, for some reason, the open magnetic flux 
tube near the surface of this pulsar is only dominated by 
the global dipolar fleld plus the multipolar magnetic fleld 
with I — 20 and m = 20, an assumption which agrees with 
the results of Mitra et al. (1999). 

The calculations according to the same scheme as in 
the previous case were done for the field lines which, at the 
initial height r. = 1.5Rs (where ^20,20 ^ Bd), constitute 
the outermost 0.1 fraction of the open field lines tube. In 
other words, the boundary conditions were taken within the 
following range: 



0.961; < e{r.) < 
< (t>{r,) < 2tt. 



(50) 



Let us notice that, just like in the previous case, angular 
spacings between the initial points within these ranges at 
the distance were taken equal, both in azimuth (f) and in 
polar angle 6. 

The magnetic fleld structure of such a configuration 
near the stellar surface is presented in Fig. |^. Similarly to 
the previous case, we observe that some field lines deviate 
outwards, that is, towards the 'appropriate' local multipolar 
poles, whereas other field lines deviate inwards. For clear- 
ness, we again do not represent in this plot (and the other 
figures below) the inner open field lines, which are less af- 
fected by the multipolar field, and hence follow the dipolar 
geometry. 

Fig. ^ represents the peripheral part of the modified po- 
lar cap of this pulsar. It is clearly seen that the outer margin 
of the polar cap appears to be fragmented into 20 symmet- 
rically arranged sub-regions. One such particular sub-region 
is plotted in detail in Fig. The azimuthal variation of 
the total magnetic field strength measured in the units of, 
first, the dipolar field, and then, the critical field, as well as 
the variation of the curvature radius of the last open field 
lines within this sub-region are displayed in Figs. pT|(a), (b) 
and (c), respectively. The variation of Pa, Bs/Bq and ps with 
the polar angle 9 towards an 'appropriate' multipolar pole 
within this sub-region is presented in Figs. 0(a), (b) and 
(c), respectively. It is seen in these figures that the condi- 
tions for the development of pair cascades are fulfilled only 
above the restricted areas of the polar cap periphery. 

At large altitudes from the stellar surface, the 'memory' 
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Figure 8. Magnetic field structure of the peripheral 0.1 fraction 
of the open flux tube near the surface of PSR B0943+10, assuming 
that the magnetic field is composed of the dipolar field and the 
multipole of order I = 20 and m = 20. 
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Figure 9. Peripheral fraction of the modified polar cap of PSR 
B0943 + 10 assuming the contribution of the multipolar compo- 
nent with I = 20, m = 20. The footprints of the open field lines are 
represented by circles. The dashed circle shows, for comparison, 
the margin of the corresponding purely dipolar polar cap. 



of the surface field structure will be preserved in tiie form of 
localized plasma columns distributed symmetrically around 
the outer margin of tlie open magnetic flux tube. Similarly 
to the previous case of Z = 4 and m = 4, the radial width of 
this thin outer ring can be estimated from Fig. |l^(c), using 
the boundary conditions it will constitute about 1/20 

fraction of the total radius of the open field lines tube. This 
'hollow cone' will also be fragmented azimuthally, and the 
transverse dimension of an individual plasma column can 
be estimated from Fig. |l^(c) as being about 1 /40 fraction of 
the circumference. These plasma columns will be separated 
by low-density 'slits' of approximately the same transverse 
angular dimension. 

Obviously, such a structure of the open flux tube ex- 
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Figure 10. A fraction of the modified polar cap (see Fig. ji]), 
representing one of the localized sub-regions where the curva- 
ture radius of the open field lines varies from 'favourable' to 'un- 
favourable' values for the pair creation process. 



hibits a striking similarity with the 20-fold pattern of emis- 
sion sub-beams discovered by Deshpande & Rankin (1999, 
2001) using their 'cartographic transformation' technique. 
They found an identical configuration at two different fre- 
quencies, which suggests a filamentary plasma distribution 
rather than an emission process, as a possible origin for such 
a pattern. Deshpande & Rankin (2001) conclude that such 
plasma columns should have feet within a certain ring on 
the polar cap. In view of our multipolar model, such a ring 
could be associated with the outer margin of the modified 
polar cap, represented in Fig. ^ 

According to Deshpande & Rankin (1999), 20 emission 
sub-beams undergo rigid counter-clockwise rotation around 
the magnetic axis in a total time of 37P = 41 s, maintain- 
ing their number and spacing despite perturbations tending 
both to bifurcate a given beam and to merge adjacent ones. 
Although at this stage it is only possible to speculate about 
several possible reasons for this effect, we suggest in the next 
section an explanation of such a rotation, in terms of the 
character of the perturbations in the azimuthal direction. 

In another way, such a phenomenon could actually re- 
fiect some stroboscopic-like effect. Indeed, pair creation in 
one of the sub-regions in Fig. ^ results in a roughly Gaus- 
sian distribution of the pair plasma density across the local 
magnetic field. Such a distribution function would screen the 
electric field, at least partially, in the adjacent sub-regions, 
thus temporarily ceasing or, at least, largely reducing the 
pair production process there. The situation would revert af- 
ter pair plasma cascade will exhaust in this sub-region, when 
the accelerating electric field is screened out. Then, two ad- 
jacent sub-regions will 'turn on', reducing pair creation in 
their adjacent sub-regions. Such consecutive 'switches', in 
combination with the E X B-drift and the polar cap heat- 
ing, may result in the above-mentioned stroboscopic effect. 

It should be noticed that, in order to obtain the pat- 
terns presented in Figs. ^ — one has to choose the free 
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Figure 11. Dependence of /3s, Bs/Bq and ps on the azimuthal Figure 12. Dependence of /3s, Bg/Bq and ps on the polar angle 0, 

angle if>, at the feet of the last open field lines of PSR B0943 + 10, in the direction of the local multipolar pole, for the same pulsar, 
assuming contribution of the multipole I = 20, m = 20. 



parameter f (introduced in equation ^) in our calculations 
so that the magnetic field strength in the close vicinity of 
the multipolar poles themselves appears to exceed signifi- 
cantly the critical field Bq. However, these multipolar poles 
are located at lower latitudes of the stellar surface, close to 
the neutron star's equator. Therefore, very high magnetic 
field is achieved well within the domain of the closed field 
lines where plasma creation processes certainly do not oc- 
cur. Besides, there are no serious theoretical limitations to 
the strength of the magnetic field in nature (Usov V. V., 
private communication) . 



4 A MODEL FOR RADIO EMISSION FROM 
AN ENSEMBLE OF FINITE SUB-BEAMS 



4.1 Features of the model 



As demonstrated in Section 3.4, the particular magnetic 



field topology which results from the superposition of the 
global dipolar magnetic field and a high-order multipolar 
field with m 7^ 0, both with high strength, suggests that the 
pair plasma created at the bases of magnetic fiux tubes is 
guided within isolated thin filaments, located at the outer 



margin of the bundle of open magnetic field lines. According 
to Sturrock (1971), the density of the pair plasma within the 
flux tubes is proportional to the Goldreich & Julian (1969) 
density, namely Up = nn^ j , the k parameter being in the 
domain 10"^ < k < 10^. Outside, the plasma is supposed to 
have a much smaller density, of the order of the Goldreich 
& Julian density, Wqj . Thus, the whole system of dense fila- 
ments can be described as an ensemble of relativistic finite 
beams immersed in external low-density media. These ex- 
ternal media in the pulsar magnetosphere are, on the one 
hand, the charged static medium lying on closed field lines, 
and on the other hand, the primary beam formed of ultra- 
relativistic particles of one charge fiowing inside the hollow 
cone and between the dense beams. 

The simplest case concerns one flowing annular beam 
of relativistic pair plasma surrounded by external plasmas. 
The dispersion relation giving the characteristics of pertur- 
bations able to propagate in such a system were obtained 
matching the solutions for the electromagnetic fleld compo- 
nents at the different beam-plasma interfaces (Asseo 1995). 
The case of the e nser nble of twenty isolated sub-beams dis- 
cussed in Section 3.5 is more complex: the beams are sup- 
posed to be regularly distributed over a ring that delimits 
the emission region; they have different radial and azimuthal 
extents and are separated over the ring by distances of the 
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order of their azimuthal extent. The radial variation of the 
density operates in two different ways: either there is an 
abrupt change of the density at the edges of the ring, or 
there is none. The azimuthal changes of the density from 
one beam to the next one are alike. The geometrical charac- 
teristics of the relativistic beams are evaluated taking into 
account the fact that the observed ensemble of twenty beams 
have their origin in the gap region, very close to the stellar 
surface, where the complex multipolar field determines the 
pair creation process. As concluded in the previous section 
(see Figs. ^2|(c) and 0(c) for details), the radial angular 
width of an individual sub-beam constitutes a fraction equal 
to about 1 /20 of the total radius of the open field lines tube 
at an altitude (where B'™ ^ B'^), whereas its azimuthal 
angular width represents a fraction equal to about 1/40 of 
the circumference. Therefore, at some distance (see Sec- 



(53) 



tion 



3.1), the radial and azimuthal extents of individual sub- 



beams can be estimated as, respectively. 



wg — 



and 



20 



7.3 X lO^P" 



3/2 



23n;^^2.3xlO^P 
40 



(51) 



(52) 



where equation ( p6| ) has been used in order to evaluate 6-. 
For example, at a distance r. = 2_Rs, far above the gap re- 
gion, the numerical estimates for the individual sub-beam 
width yield wg ~ 2 x 10^ cm and ~ 6.5 x 10"^ cm, 
respectively, for a pulsar with period P = 1 s. The ratio 
of the azimuthal and radial sub-beam widths is constant, 
namely w^/wg ~ tt whatever the distance, but the typical 
wavelengths associated with observed radio frequencies are 
smaller than these characteristic widths. Thus, one should 
test whether it is sufficient to consider two-dimensional per- 
turbations of the whole system of beams immersed in ex- 
terior media, or if it is necessary to analyse the fate of 3- 
dimensional perturbations. This can be checked as in Asseo, 
Pellat & Sol (1983), considering the strongest instabilities 
which may develop in a thin beam surrounded by difi'erent 
plasmas. 

In the previous sections we used the spherical system of 
coordinates {r,9,(j>), centered at the stellar centre. Instead, 
in this section it is convenient to use a cylindrical system of 
coordinates {q,^p,z), adequate for a local description. Here 
the radial coordinate g is measured along the local curvature 
radius of the magnetic field line; (p is the azimuthal angle 
which corresponds to the angular variation along the circle 
with the radius p, this circle being tangent to the magnetic 
field line passing through the point (r, 0) in the plane (j) = (po 
of the spherical geometry; z is the coordinate perpendicular 
to the plane 4> = 4>o- 

Assuming 3-dimensional perturbations in the cylindri- 
cal system depending on [q, 93, z) and characterized by wave- 
numbers m and kz , the strongest instabilities are obtained at 
the resonant frequency u ~ mc/ p and for transverse wave- 
numbers in the domain Akz < yplrr?^^ p~^ . In this case the 
analytical solutions obtained for the electromagnetic fields 
characterizing the perturbation in terms of Bessel and Han- 
kel functions, result in an unmodified growth rate of the 
instabilities. Consequently, they are valid for beams with a 
spatial extent 



where we used the fact that the spatial extent of a sub- 
beam Az = 2-K / Akz is nothing else but its azimuthal width 
ui0 (equation M). 



4.2 The case of millisecond and fast pulsars 

Numerical estimates show that in the case of millisecond 
and fast pulsars the constraint ( p3[ ) is easily fulfilled just 
above the gap region and beyond, considering the domain 
of observed radio frequencies v = 2-k/uj (see Figs, p^a) and 
(b), respectively). This indicates that in such a case, finite 
transverse dimensions do not modify the results obtained in 
cylindrical geometry using two-dimensional perturbations to 
characterize radio emission features generated above the gap 
region within an annular hollow conical beam, or within an 
ensemble of discrete sub-beams similar to those observed in 
PSR B0943 + 10. 

In such conditions, one may consider that locally the 
propagating wave only sees the radial variations of the den- 
sity, so that it is possible to treat the problem in a poloidal 
plane. While doing so, we ignore the density variations in 
the direction corresponding to the azimuthal (/)-direction of 
the spherical geometry, that is to say, in the transverse z- 
direction of the cylindrical geometry, and treat each beam 
of the ensemble of beams as isolated and immersed in 'in- 
finite' exterior media. This allows us to extend the results 
which concern the characteristics of the perturbations able 
to propagate in the system formed by relativistic beam and 
plasma flows which fill a hollow conical region and are im- 
mersed in exterior media, to the more complex situation of 
the ensemble of beams. 

The dispersion relation specific to the case of a finite 
beam immersed in external media was obtained by match- 
ing the solutions for the electromagnetic field components 
at the different beam-plasma interfaces (Asseo 1995). It can 
be expressed in terms of the dielectric constant for the rel- 
ativistic beam plus relativistic plasma flows, W , and of the 
dielectric constants of the exterior media, Wi and We- It is 
obtained close to the resonant frequency. 



(54) 



where mf2o ~ mc/ p is the frequency associated with the 
circulation of relativistic particles along magnetic flux tubes 
and Y = Sui/ui^ < 1 is the relative growth rate. Assuming 
variations with the radial distance together with a wave- 
like behaviour proportional to exp [—i{uit — rrKj))] for linear 
perturbations of the whole system, the dispersion relation 
becomes: 



tanh • 



2mW^^^ 



-Y) 



3/2 



3/2' 



w 



1/2 



(55) 



Let us notice that the radial width in the cylindrical sys- 
tem of coordinates, Wg is approximately equal to the angular 
width Wg in the spherical geometry, so that the above disper- 
sion relation is correct in a cylindrical geometry. Such a dis- 
persion relation for finite beam and plasma flows differs from 
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Figure 13. The condition (^) at different distances from tfie 
stellar centre. The panels (a), (b) and (c) represent the results of 
calculations for the sample pulsars with periods equal to 1.6 X 
10-3 s, 0.1 s and 1.0977 s (PSR B0943 + 10), respectively. In each 
of the panels, solid curves correspond to the observed frequency 
u = 34 MHz, dashed curves to v = 111.5 MHz, dashed-dotted 
curves to v = 430 MHz and the dotted curves to v = 1.4 GHz. 



the dispersion relation that characterizes the two-stream in- 
stabihty obtained for infinite and homogeneous beam and 
plasmas described in a straight geometry, namely W = G. 
It is modified due to the presence of external media with 
dielectric constants Wi and W^. It is also modified due to 
geometric properties of the flow: more precisely the curva- 
ture of the magnetic field, as the radius of curvature p of 
magnetic field lines is involved, the finite width of the beam 
We and the azimuthal wave-number of the perturbation m. 
The analysis of this modified dispersion relation in different 
limits shows that it allows the description of both the ra- 
diative (see Goldreich & Keeley (1971) who first described 
it) and two-stream instabilities, well-known in the context 
of pulsar physics, but also that there is an unexpected insta- 
bility: the 'finite beam instability', which even develops for a 
one-component beam of relativistic mono-energetic particles 
bounded by external media. 

The characteristics of the 'finite beam instability' were 
obtained in the most general conditions for a beam bounded 
by identical, or non-identical, external media from the 



above dispersion relation (Asseo 1995). The case of a beam 
bounded by identical external media with Wi = M^e, is of 
interest here for an ensemble of beams delineated by the su- 
perposition of a strong multipole magnetic field plus a dipole 
field. In the 'thin' beam case with Y ^ we / p, the dispersion 
relation simplifies to: 



\Wi) 



1/2 



tanh 



£^2 



1/2^1/2 



(56) 



The dependence on the presence of external plasmas appears 
through Wi. The dependence on the geometric properties of 
the flow appears through the parameter: 



m we 



(57) 



For the magnetic configuration which includes both 
dipolar and multipolar magnetic field components, the value 
of the radius of curvature of magnetic field lines varies with 
the distance from the stellar surface. Within the gap re- 
gion from which the beams are issued, the radius of cur- 
vature of magnetic field lines is of the order, or even less, 
than the radius of the star, p < Rs- Beyond, in the re- 
gion where the dipolar magnetic field is dominant, the ra- 
dius of curvature is approximately p ~ 10^ -Rs, (indeed, 
p ^ (Ar/Se) ^ 10** cm at the distance r = 2Rs). There- 
fore, the azimuthal wavenumber which corresponds to the 
resonant frequency at which the instabilities are triggered 
off, fixed by the relation m ~ pu>j^/c, slightly varies in the 
radio emission domain which starts above the gap region 
(see Asseo 1995). At such resonant frequencies. 



£1. 

m We 



c 1 

Lj We 



(58) 



only depends on the radial extent of the involved annular 
beam or sub-beams. Numerical estimates with usual param- 
eters for a pulsar with P = 1 s show that in the domain of 
observed radio frequencies lo^ < 1.6 x 10® Hz and for the 
radius of curvature of magnetic field lines in the domain 
p < 10* cm, e is smaller than unity. 

Above the gap region, but relatively close to it, the ratio 
of the ring width to the radius of curvature of magnetic field 
lines is very small, we j p < 2x 10"^, for values of the radius of 
curvature in the domain p < 10* cm. Thus, we are in the thin 
beam case, the dispersion relation (|5^) is valid in its limiting 
form (equation B6|), and the approximation tanh{. . .} « 1 is 
adequate, the parameter e being much smaller than unity. 
In this case, the dispersion relation (^^ simplifies to W — 
— Wi. Such a dispersion relation is adequate to describe a 
limited region of emission whatever its altitude. 

For a one-component beam of relativistic particles with 
Lorentz factors 7p, the dielectric constant W is defined as: 



W = 1 



y2 ' 



(59) 



and the dielectric constant Wi of the interior plasma sup- 
posed to be at rest is defined as. 



m = 1 
Here, 



Vtir)il-2Y). 



3/2 r> ' 

7p' milo 



(60) 



(61) 
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whereas ^tJp{r) — [47rnp(r)e^/m] ^''^ and ujpi{r) = 

47rnpi(r)c^/m] ^''^j the plasma frequencies in the flowing 
beam and in the exterior media, respectively, vary along 
magnetic flux tubes with the corresponding densities. Thus, 
■p(r) and 'Pi(r) depend on the densities in the interior and 
exterior flux tubes, respectively. Assuming that, as men- 
tioned above, the density within the magnetic flux tube is 
proportional to the Goldreich & Julian (1969) density, 

np{r) ^ KU^jir), (62) 

whereas the density in the exterior media npi(r) is equal to 
the Goldreich & Julian density, 

-oAr) - ^ - 7 X 10-§ [cm]-, (63) 

we have, 

W»L5xlO«^4(ff-j4^. («) 

V.ir) » 1.5 X U'^I^^J^. («) 

Let us note that Wi varies according to the value of 
■Pi(r) in the emission region, namely: 

(i) Wi ~ 2Y at the transition zone located at the distance 
r « r^.^ where 'Pi(r) ~ 1. 

(ii) Wi ~ ~'P'i{i') close to the surface of the star where 
'Pi{r) S> 1, at distances r < r^^. 

(iii) Wi ~ 1 beyond the transition zone where Viij-) <^ 1, 
at distances r > r^^^. 

Including a multipolar magnetic field superimposed 
over the dipolar field and assuming that at the surface of the 
star the intensity of the multipolar magnetic field is propor- 
tional to the strength of the dipole field^Bg"' = (/3s - l)B'^ 
(where (3a has been defined in equation Bll) , we have for the 
total field, B{r) = Bi{RJrf + (/J^ - l)Bf(i?s/r-)'+^ There- 
fore, the distance r^.^, at which Viir) « 1, is not significantly 
different from the distance r^.^^ obtained in the case where 
there are no multipolar field components in addition to the 
dipolar field, as the relative contribution of the multipolar 
magnetic field, {(3^ — l)(-Rs/»')'~^, is very small. In effect, r\^^^ 
is simply determined from the relation 

'P^ir,,J=xfir,^J = l, (66) 



where 




Tjjj^ increases with the strength of the total magnetic field 
but decreases as the pulsar period and the frequency as- 
sociated with the circulation of relativistic beam particles 
increase. 

In the extreme case of a millisecond pulsar with pa- 
rameters P = 0.0016 s, Bt = 10** G, mfio = 10^ Hz, we 
obtain r^.^ = 2.4 x 10^ cm. For a 'normal' pulsar with 
parameters P = 0.1 s, B = 10^^ G, mfio = 10^ Hz, we 
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obtain r^.^ ~ 6 x 10 cm. For a 'slower' pulsar with pa- 
rameters P = 1 s, B = 10^^ G, mf^o = 10^ Hz, we obtain 
''lim ~ 3 x 10® cm. Therefore, whatever the type of pulsar. 

The distance r^.^^ is important in delineating the differ- 
ent regions of instability: below r^.^ the 'finite beam insta- 
bility' is dominant whereas beyond r^.^ either the radiative 
or the two-stream instabilities prevail. As a matter of fact, 
the behaviour and growth rate of the 'finite beam instability' 
is different according to the distance of the considered emis- 
sion region to the stellar surface. The characteristics of the 
unstable waves excited trough the 'finite beam instability' 
- frequency, polarisation, available electromagnetic energy 
and growth rate - are in agreement with pulsar radio ob- 
servations (Asseo 1995). Thus, the 'finite beam instability' 
will initiate the radio emission process in the ensemble of 
beams with their origin in the gap region, relatively close to 
the surface of the star. Moreover, such an instability strongly 
depends on the ratio of densities at the basis of the magnetic 
tube within and outside the limited beam. As the density 
of the created pair plasma may vary, an observed apparent 
motion of the beams around the magnetic axis could be re- 
lated to variations of the densities within and outside the 
successive tubes. 

4.3 The case of slow pulsars and PSR B0943 + 10 

Numerical estimates concerning the above constraint on the 
spatial extension of the beams, required to ignore azimuthal 
variations of the perturbations (equation B3) , show that in 
the case of PSR B0943 -I- 10 this constraint is not uniformly 
satisfied in the whole region above the gap (see Fig. p^c)). 
The distance r^^ from the center of the star beyond which 
the constraint is satisfied varies with the frequency of ob- 
servation: as clear from Fig. |l^(c), larger distances corre- 
spond to lower frequencies. Indeed, r,^ ~ 2.3 x 10*' cm 
for a frequency of observation u = ui/2n ~ 430 MHz, 
^ 4.6x lO'^ cm for ^ 111.5 MHz and r^^ ^ 8.7x lO'' cm 
for u^34 MHz. 

Consequently, between the gap region and the loca- 
tion r — r,Q, it is necessary to take into account the 3- 
dimensional character of the perturbations able to propagate 
in the system of discrete sub-beams: the density variations 
in the direction corresponding to the azimuthal 0- direct ion 
of the spherical geometry, that is to say, in the transverse 
z-direction of the cylindrical geometry, cannot be further ig- 
nored and the ensemble of sub-beams has to be considered as 
a whole system immersed in 'infinite' exterior media. Thus, 
an additional systematic azimuthal variation of the pertur- 
bations below will exist and modify the character of per- 
turbations beyond r^^ , in the domain where 2-dimensional 
perturbations are dominant. In fact a coupling of the differ- 
ent perturbed physical quantities is involved as soon as one 
introduces an azimuthal dependence of the perturbations at 
relatively small distances from the surface of the star. A 
thorough treatment including the dependence on the coor- 
dinates (r, <j), z) of the specified cylindrical geometry of the 
perturbations is possible as in Asseo et al. (1983), writing 
the perturbations as of the type: 

g{z, r) exp [—ik^z] exp [—[{ut — m(/))]. (69) 
Obviously, the character of the perturbations in the region 



18 Asseo & Khechinashvili 



located beyond r^^ has to approach the one described above 
for fast pulsars, whereas in the region below r^j^ it is modi- 
fied by the dependence on the variable z and wave-number 
kz- Although we here restrict to a qualitative explanation, 
and leave exact calculations concerning the character of 3- 
dimensional perturbations for future work, this suggests that 
for PSR B0943 + 10, the observed 'rotation', and/or 'drift- 
ing', of the emitting columns may be a trace of the initial 
waves formed in the region below r^^ and of their transition 
towards the above region. 

Most of the pulsars which exhibit beautiful drifting sub- 
pulses - PSRs 0809 + 74, 0320 -|- 39, 0820 -|- 02, 0818 - 13, 
1540 — 06 and 0943 + 10 show periodic drifts (as observed 
and/or classified by Rankin (1986), sec table 2 therein). We 
notice that all the quoted pulsars have sufBciently long pe- 
riods to be modelled in analogy of PSR B0943 + 10. This 
suggests that 'drifting' effects, already known as intrinsically 
conal phenomena and associated with the circulation of dis- 
turbances within the annular polar cap region, could also 
result from similar perturbations moving within a hollow 
conical region fragmented into discrete sub-beams. Then, 
azimuthal widths of the sub-beams will be reflected in the 
character of unstable excited perturbations. 

5 SUMMARY 

In this paper we study the complex magnetic topology in 
the vicinity of neutron stars. We restrict to the case of the 
aligned rotator, assuming parallel, or anti-parallel, orienta- 
tion of the rotation and magnetic axes. We assume that the 
total magnetic field results from the superposition of the 
global dipolar field and of star-centered multipolar fields of 
comparable, or even of higher strength. We obtain the fol- 
lowing results: 

(i) For configurations involving axially symmetric and 
uniform multipolar magnetic fields it is not possible to ob- 
tain small curvature radii of the magnetic field lines, even 
for multipoles of higher order and strength relative to the 
dipolar field. Let us recall that according to pulsar polar gap 
models, curvature radii of the order of the stellar radius are 
crucial for an efficient pair creation process. 

(ii) An interesting consequence of such a magnetic topol- 
ogy is the existence of 'neutral points' and 'neutral lines' 
above one of the pulsar polar caps. We plan to study impli- 
cations of such peculiar features of the magnetic field geom- 
etry in future works. 

(iii) For configurations involving axially symmetric and 
non-uniform multipolar magnetic fields, the magnetic field 
structure just above the polar cap is modified in such a way 
that the pair plasma flow at high altitudes appears as modu- 
lated in separate filaments regularly distributed around the 
outer margin of the open flux tube. Such fllaments have their 
foot in the modified polar cap. They are thin and have finite 
radial and azimuthal extents. Therefore, the non-uniform 
pattern of the distribution of curvature radii over the po- 
lar cap is refiected in the inhomogeneous distribution of the 
pair plasma density within the cross-section of the open flux 
tube at high altitudes. 

(iv) For the particular superposition of a dipolar field and 
a multipolar field of order I = 20 and m = 20, the mag- 
netic topology causes fragmentation of the pair plasma flow 



into twenty isolated thin filaments. This resembles the set of 
twenty sub-beams observed by Deshpande & Rankin (1999, 
2001) in PSR B0943 + 10. 

(v) Emission in such a system of isolated thin beams, 
fiowing along the curved magnetic field lines, can be de- 
scribed in terms of unstable waves excited by specific in- 
stabilities. Relativistic finite beams flowing in the pulsar 
magnetosphere along a limited fraction of the hollow cone 
surrounding the magnetic axis can be treated as immersed 
in exterior media, namely the medium in the closed field 
line region and the medium within the hollow cone of the 
open magnetosphere (Asseo 1995). Three different instabil- 
ities may develop in such finite beams, depending on the 
distance to the center of the star: the 'finite beam', ra- 
diative and two-stream instabilities. As previously obtained 
by Asseo (1995), close to the surface of the star the 'finite 
beam' instability will generate waves whose properties (fre- 
quency, polarisation and available electromagnetic energy) 
are in agreement with the observed pulsar radio emission. 

(vi) We check that the conditions for such a treatment 
are relevant for the configuration of discrete sub-beams that 
may exist in millisecond and fast pulsars. They are also rel- 
evant for slower pulsars, like PSR B0943 + 10, beyond some 
distance from the gap region. On the other hand, closer to 
the gap of such a slow pulsar, the radial and azimuthal vari- 
ations of the fiowing sub-beams cannot be separated; the 
system of sub-beams has to be treated as a whole; the 3- 
dimensional character of the perturbations will introduce 
some 'drift' of the emitted waves; such a 'drift' effect will be 
transferred to high altitudes due to the necessary continuity 
between the vicinity of the gap and the farther region. In 
this way, the continuity of the fiow all along magnetic fiux 
tubes and of the perturbations carried by it, accounts for 
the 'drifting' character of the subpulses observed in pulsar 
radio profiles. 
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